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A novel model of discretized energy cascade generated by Benjamin-Feir instability is presented. 
Conditions for appearance of direct and inverse cascades are given explicitly, as well as conditions 
for stabilization of the wave system due to cascade termination. These results can be used directly 
for explanation of available results of laboratory experiments and as basic forecast scenarios for 
planned experiments, depending on the frequency of an initially excited mode and steepness of its 
amplitude. 
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Benjamin-Feir instability (BF-instability) is one of fun- 
damental principles of nonlinear water wave dynamics, 
[1]. This phenomenon is of the utmost importance for 
description of dynamics and downshifting of energy spec- 
trum among sea surface waves, formation of freak (or gi- 
ant) waves in oceans and wave breaking. BF-instability 
and its physical applications have been profoundly stud- 
ied during the last few decades, [2-5], and its main fea- 
tures can be briefly summarized as follows: 

I. initial exponential growth of the main side bands of 
the carrier wave; 

II. late asymmetry of sidebands and temporary fre- 
quency downshifting; 

III. discretized spreading of energy to higher and lower 
frequencies; 

IV. existence or absence of near recurrence Fermi- 
Pasta-Ulam (FPU) phenomenon in no breaking regime 
of wave propagation, in experiments with different pa- 
rameters of initial excitation. 

As for surface gravity waves on deep water resonant 
interactions occur at the third order, BF-instability can 
be described at early stages of the process as interaction 
of three monochromatic wave trains: carrier (wc ), upper 
(w+ — ujc + Aw) and lower (w- = lJc — Acj) side-band 
waves with Aw > which form a resonant quartet for 
one particular configuration which occurs when two of 
the waves coincide, with frequency resonance condition 



(1) 



where uji = w(ki) and are notations for dispersion 
function and wave vector correspondingly. 

In this Letter we present a novel model of BF- 
instability basing on the model of dynamical cascade gen- 
eration in a wave system with narrow initial excitation, 
first introduced in [6] . We demonstrate that in the frame 
of our model a) main features I. -IV. of BF-instability are 
naturally reproduced; and b) dependence of the cascade 
form on details of initial excitation (choice of frequency 
and wave steepness) is in accordance with available ex- 
perimental data. 



Dynamical system corresponding to an isolated quartet 
reads: 



i Ai = TAIA3A4 + (wi - uji)Ai , 

i Ai = TAIA3A4 + (W2 - UJ2)A2 , 

i A3 = T*AlAiA2 + (w3 - ^3)^3 , 
i Ai = T*A'^AiA2 + {Qi - uJi)Ai , 

^j~^j^j:t=imAi\'-'2Tn\A'^ 



(2) 



and 



where interaction coefficients Tij — Tji 
T = are responsible for the nonlinear shifts of fre- 
quency and the energy exchange within a quartet cor- 
respondingly; (wj — ujj) are Stokes-corrected frequencies 
and Aj are slowly changing amplitudes of resonant modes 
in canonical variables. Explicit form of interaction coeffi- 
cients for surface gravity waves is given in [16]. Analytical 
solution of (2) in terms of elliptic functions is found in 
[7] and is studied numerically. 

Aiming to analyze quartet dynamics qualitatively in 
the case when only one or two modes are initially ex- 
cited, one can use standard change of variables Aj = 
—iCj exp[— i(^j] , and rewrite (2) in amplitude-phase pre- 
sentation as 

dCl 
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dCl 



dt dt dt dt 

= 2 1 T I Ci C2 C3 C4 sin(arg T - 



</'12,34), (3) 



where the dynamical phase (^912.34 = (fi + (f2 — ^3 — 
corresponds to the chosen resonance conditions of the 
form 



Wl + W2 = W3 + ^ii ki + k2 



(4) 



Sys.(3) has three independent Manley-Rowe constants of 
motion 



/i3 — + C3 , /i4 — Ci + Ci , I23 — C; 



cl 



(5) 



Any linear combinations of these three are also constants 
of motion, e.g. 

I24 — C2 + cl , /1234 ~ Cl + C2 + C3 + C4 , (6) 



Ii,2 — Cl — C- 



2 ; 



-^3,4 — C3 — C4 



(7) 
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and can be used for qualitative analysis of a quartet dy- 
namics for specific initial conditions. 

Case 1: one mode is initially excited. Let initially the 
mode wi be excited, i.e. at the time moment t = we 
have 

Ci,o ^ (^2,0 — C'a.o — Cifi = Co , (8) 

where Cj^ = Cj{t = 0). It follows from (5), that the 
amplitudes C2{t), C^^t) and C4{t), being initially small, 
remains small at all time moments t > 0. Indeed, 

2co ^ hi + In — h,i = C2 + 

^ C2(t)-C3(t)^Co, (9) 
0~/l3-/l4 = C|-C| ^ C4(i)~C0, (10) 

for any t > 0, i.e. the energy transfer in this case is 
essentially suppressed. 

Case 2: two modes are initially excited. There exist 
two different types of modes' pairs which should be re- 
garded separately, [6]: one-side-pair and two-side-pair, 
referring to resonance conditions of the form (4). Ac- 
cordingly, there are two 1-pairs of modes' frequencies: 

1-pairs: (wi,W2), ((^3,^4), (11) 

and four 2-pairs of modes' frequencies: 

2-pairs: (^1,^3), (^1,^4), (w2,W3), (w4,W3). (12) 

Case 2a: 2-pair is initially excited. Let 

C'l.o — Cs.o ^ C2,o — (^4,0 = Co , (13) 

then again it follows from the form of Manley-Rowe con- 
stants of motion that C2{t) ~ C4{t) ~ cq. 

Case 2b: 1-pair is initially excited. Regard 1-pair 
{coi , u)2) with initial modes' amplitudes as follows: 

Ci,o ~ C2,0 > C3,0 ^ ^4,0 • (14) 



where 



Ai,o 


= "^'-'1,0 + -'12(-'2,0 
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A4,0 


— -'4l'-^i,o + ^4202^0 > 




W3 


= (Ai,o + A2,o - A3,o 


+ A4,o)/2 


UJ4 


= (Ai,o -1- A2,o + A3,o 


- A4,o)/2 



-?2 = i^p-^EM' (1^) 

= \T\^cl^cl^ - (riC?,o + r2C|,o) 74 , 

V = T*Ci,oC2,oexp[i(Ai,o + A2,o)t], 
Ti = ^Tn + T12 + Ti3 + Ti4 , 

72 — T12 + 2^22 + 723 + 724 • 

Accordingly, the evolution of the amplitudes ^3 and A4 
is defined by the sign of the increment :^^2 given by (16). 

If > 0, (16) predicts exponential grow of ampli- 
tudes As and A4. In this case energy goes from the ini- 
tially excited 1-pair (cui , 1^2) to the second 1-pair (073 , W4) 
with characteristic time ~ l/z^i2- Similarly, if the 1-pair 
(cjs ,0^4) is initially excited, energy can effectively go to 
the 1-pair (wi ,0^2) , if ^'34 > 0. 

However, (16) does not guarantee that 1/^2 is posi- 
tive. For instance, if Ci,o ^ C2,o or Ci,o -C C2.0, ^12 is 
negative for any interaction coefficient while in the case 
C\fi = 6*2,0 the sign of j/^2 depends on the relations be- 
tween \T\ and (71 + 7^)/2: 

^^?2 = [m'-(ri + r2)V4]cf,o, (17) 

If < 0, the solution (15) yields pure oscillatory 
behavior of the amplitudes A3 , A4 with frequencies a73 ± 
\ui2\ and u)4 ± |z^i2| correspondingly. 



There exists no restriction on modes' growth originating 
from the Manley-Rowe constants of motion, and result- 
ing evolution depends on the details of the initial energy 
distribution within a quartet. 

During initial evolution, during which inequalities (14) 
still hold, one can neglect the effect of modes ^3,0^4 on 
modes a;i,W2 and Eqs. (2) can be solved explicitly: 

!Ai(i) =Ci,oexp(iAi,ot), 
^2(t) =C2,oexp(iA2,of), ^^g^ 
Aait) = C3,oCxp [{ius + i^i2)t] , 
Al{t) = 6*4,0 exp [{-1104 + '^i2)t] , 



3. GENERATION OF A CASCADE 

The general model of dynamical cascade generation in 
3- and 4-wave systems with narrow initial excitation is 
sketched in [6]. In this Letter we work out the details 
and apply it for the description of BF-instability. 

At the initial step, n = 0, frequency resonance condi- 
tions have the same form (1) for both direct and inverse 
cascade while the only form of a quartet in which one- 
mode excitation yields generation of resonant interactions 
reads 

Wl -I- W2 = 2W3 (18) 
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and only in the case when the mode with frequency W3 
is excited. This occurs due to Hasselmann's criterion 
of instability for 4- wave systems, [8]. For all other con- 
figurations and initial excitations, any single mode in a 
quartet is neutrally stable. The most effective resonance 
takes place if (18) is satisfied exactly, i.e. Aw should be 
the same for both a;+ and . 

Hasselmann's criterion can be applied at each further 
step yielding the general form (18) at each step of a cas- 
cade. To simplify further presentation we introduce no- 
tation U)±n = Wc ± nAoj. 

At the step 1, a couple of new modes coi and oj-i is 
generated. From pure kinematical considerations, at the 
next steps of a cascade all possible quartets of the form 

Wc ± =F t^-i = '^c ± 2Aa; (19) 

are possible. However, different dynamical properties of 
1-pair and 2-pair in a quartet define two possible com- 
binations only, yielding exact frequency resonance condi- 
tion for quartets of the form (18) with excited 1-pair of 
modes: 

Wc -|- Wc + 2A(x) = 2{ujc + Aw) 

Wc + a;2 = 2wi, (20) 
cOc + oJc — 2Aui = 2{ujc — Aw) 

Wc -I- w_2 = 2a;_i. (21) 

Accordingly, the beginning of direct and inverse cascades 
is given by (20) and (21). As Wc-|-w±2 = 2w±i, frequency 
shift |wc — w±2| = 2 Aw occurs (its sign is opposite for 
direct and inverse cascade). 
Similarly, at the step n we have 

Wc + Wc + 2nAw = 2(wc + riAw) 

^ w„+i -I- w„_i = 2w„ (22) 
ijOc + oJc — 2nAw = 2(wc — nAw) 

^ ^-{n+l) + ^-(n-1) = 2w_„, (23) 

and the frequency shift is again 2Aw, i.e. its magnitude is 

the same at each step of both direct and inverse cascade. 
The complete system describing both direct and inverse 
cascade takes now the following form: 

W-i- + w_ = 2wc 
Wc -h W±2 = 2w±i 
< w±i -I- w±3 = 2w±2 (24) 

,'^±(n+l) +'^±(n-l) = 2W±„. 

Direction of cascade depends on the sign chosen in the 
lower index w±„. Indeed, the choice of 

uj—n generates a 

sequence of frequencies Wc > w_i > ... > w_„, i.e. inverse 
cascade, while the choice of w+„ yields direct cascade 

Wc < Wi < ... < W„. 



4. DYNAMICS OF CASCADE 

Exponential growth of the two main side-bands Wi, w_i 
at the initial step n = is caused by the resonance (18) 
at the expense of the initially excited mode W3 = Wc and 
has non-dimensional increment (in physical variables, i.e. 

< 7 = £^(2 - <t2)V2/2 < 1, (25) 

where a = Aw/eujc is the ratio of the relative resonance 
frequcincy bandwidth of the wave train to the initially 
excited wave steepness e = Cckc- 

Maximal growth rate is the same for sub- and super- 
harmonics of the excited wave: 

Imax = eV2 = C%y2 (26) 

At the next step n = 1 direct and inverse cascade will be 
be initialized by two already exciting modes Wi and w_i 
correspondingly. 

However, the important observation is that at the first 
step of a cascade the maximum of the increments for 
direct /dir,+ and inverse /inv,- cascades differ: 

/dir,+ ~ Clkl/2 7i„v,- ~ CHl/2. (27) 

This means that in the regime of symmetrical growth of 
the main side-bands C+ ~ C_, the asymmetrical behav- 
ior of the growing modes will be observed since 

k+> k- => /dir,+ > /inv -• (28) 

Similar considerations allow us to conclude that at the 
n-th step of the cascade again 

/dir,+n ^ 7inv,— n* (29) 

The last equation shows that at each step, increment of 
instability for direct cascade is bigger than increment of 
instability for inverse cascade. 

This clearly explains the frequency downshift phe- 
nomenon - the main higher side-band component spreads 
energy to higher frequency modes 

Wc < Wi < ... < W„ 

by direct cascade mechanism much faster than the main 
lower side-band by the inverse cascade with correspond- 
ing sequence of frequencies 

Wc > W_i > ... > W_„. 

Accordingly, main lower side-band mode finally will be 
dominant which is manifestation of frequency downshift 
phenomenon. 

This scenario can be clearly observed e.g. in [5] where 

laboratory observations of wave group evolution, includ- 
ing breaking effects are presented. In particular, it is 
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stated there (p. 223) that "the initial distribution of en- 
ergy (...) is altered by the transfer of energy to free waves 
u)o±Su), first noticeably to n = +2, then n = —2, +3, etc. 
These relatively fast transfers seems to be a consequence 
of detuned resonances of which first is 



I (wo) T ('^o — Suj) ± {ujQ + Suj) = cjo ± 26oj, 
^{ko)T {ko - 5k) ± (fco + Sk) = ka± 25k + Afc, 



(30) 



in which Afc is a small detuning factor." As in our no- 
tations Wo = Wc and = Aw, the frequency resonance 
condition above is just a more compact form of two fre- 
quency resonance conditions (20), (21) for quartets of the 
form (18) with excited 1-pair of modes. 

Experiments presented in [5] have been performed in 
a large wave tank (50 x 4.2 x 2.1 m), for the range of 
carrier wave initial steepness 0.12 < e < 0.28 and wave- 
lengths 1 to 4 m. Both direct and inverse discretized en- 
ergy spectra have been clearly observed (see e.g. Fig. 20 
from [5] ) , breaking and recurrence phenomena have been 
reproduced. 

The authors concluded that " a serious challenge is im- 
posed by these results to any predictive method of ocean 
wave evolution in which downshifting depends solely 
on slow, high-order, resonant wave interactions" which 
"might very well dominate the downshifting of energetic 
waves in the ocean, too" ([5], p. 226). 

We conclude here that Tulin and Waseda in [5] came 
very close to the conclusive solution of this problem while 
describing the first few steps of dynamical cascade. Only 
one final step - realization that corresponding detuned 
resonances build clusters of the form (24) - has not been 
done. 



5. TERMINATION OF CASCADE 

Energy cascading will be terminated as soon as the 
BF-criterion of instability (25) is violated, i.e. if / w 
(stabilization of the wave system due to transition to 
linear regime) and if / > 1 (substantial increasing of the 
cascading modes' steepness and transition to breaking 
regime). In the last case 



cr^ = (Aw/wc/e)^ > 2. 



(31) 



For the n-th step of cascade, the magnitude of can be 
estimated as follows: 



(w„e)2 (wc + nAujfCl (fce + nAfc)^ 



(Aw)^ 



>2. (32) 



Wc can obtain upper estimate of the energy spreading 
at the n-th step of the cascade in the resonance of form 
(18) by its maximal value, when the excited carrier wave 
energy will be totally and equally distributed between 
two resonance side-band modes: 



This means that the general upper estimate for the n-th 
step of cascade has the form 

Cl < GlIT. (33) 

After substituting (33) to the condition (32), this yields 



2"(Aw)2 



(Aw)^ 

(a;„e)2 (w^ + nAw)2C2 (fc^ -h nAfc)2 



>2. (34) 



The latter inequality shows clearly that the stability 
criterion will be reached after a finite number of steps 
for the direct energy cascade; correspondingly, energy 
spreading to higher frequencies will be terminated due 
to the growth of nonlinearity and consequent breaking 
effects which is in accordance with laboratory results [5] . 

For the inverse cascade we can estimate the stabiliza- 
tion condition in a similar manner as 



(Aw) 



2" (Aw) 



(w_„£)2 (wc - nAw)2C^(fc^ - nAfc)2 



>2, 



(35) 

i.e. the inverse cascade will be terminated even faster 
than the direct cascade. 

These estimates have been made in accordance with 
classical results (25), (26) and are valid only for small 
enough initial steepness of the excited wave e ~ 0.1. 

For moderate and high initial wave train steepness 
e ~ 0.1 0.4, improved results for the growth of incre- 
ments and a criterion of instability have been obtained 
by Dysthe, [9] . In this case, instability violates if 



9 2Aw\ 
(T^ > 2(1 ) 

V Wc / 



(36) 



Cl+, < C'J2. 



This obviously means that termination of energy cas- 
cades can happen even earlier for a higher initial steep- 
ness of carrier wave. 

As it is shown in [9], maximal growth rate is modified 
by the following way (see Eq.(3.10), p. Ill, notations as 
in the cited paper): 

7„ = i(l-2Ao)^2 

where = ak is steepness of carrier wave. Accordingly, 
in our notations 

Ima. = ^£'(1 - 2e) = ^C^il - 2Ccfcc). (37) 

Expressions for /dir,-i-n and /inv,-n at the first cascade's 
step read 

-fdir,+ = ^C+fc+(l - 2(7+fc+), fc+ = fce + Afc, (38) 
/dir,- = ^C2fc2(i-2C_fc_), fc_=fcc-Afc, (39) 
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and in the regime of symmetrical growth of the main 
side-bands C+ ~ C_ ~ C we have 

Idir,+ - Idir,- ~ 

^C^[kl{l - 2Ck+) -kl{l- 2Ck_)] 

= lc'[{kl-k'_)-2C{kl-e_)] 

= ^C^{k+ - k-)[k+ + k^- 2C{k\ + fc+fc_ + k^_)] 

= 2C^Ak[kc - C(3fc2 + (Afc)2)].(40) 

As C, k, Ak arc positive, it follows from (40) that idir,+ > 
Idir,- only if the following condition holds 

kc - C{3kl + {Akf) > 0, (41) 

which is satisfied in the finite range of wavelengths kc- 
Similar considerations show that at n-th cascade's step, 
-^dir,-i-n > -'^dir.-n IS also Satisfied in the finite range of 
kc such that kc - C{2>kl + (nAkf) > 0. For all other 
kc, maximal increment for inverse cascade will be larger 
than for direct cascade, /inv,-n > -'^dir.+n- 

Thus, by considering two forms of instability increment 
given by (26) and (37) for the energy cascade concept, we 
can make the conclusion that for essentially high initial 
steepness e > 0.25 the inverse dynamical energy cascade 
to lower frequencies may dominate. In general, both di- 
rections of cascade may be significant, unlike the case of 
a small initial steepness of the excited wave. 

6. CONCLUSIONS 

• The main features I.-IV. formulated in the Introduc- 
tion have been previously studied mainly numerically, 
while our model gives a clear explanation of the physi- 
cal origin of the observed phenomena. 

• Moreover, several fundamental facets of BF-instability 
that have not been investigated theoretically until now - 
cascade direction, its dependence on the initial conditions 
and finiteness of the number of cascading modes - are 
adequately and constrTictivcly described by our model: 

•• Direction of cascade depends on the initial wave 
steepness of the wave train; for small enough steepness 
e ^ 0.1 direct energy cascade prevails, while for high 
enough steepness e > 0.25 inverse cascade and energy 
spreading to lower frequency modes may be comparable 
with direct cascade rates. 

•• An increase of initial steepness may reduce the num- 
ber of steps for energy cascades in both directions. 

•• Frequency downshifting in the wave system is 
caused by the direct energy cascade and spreading of en- 
ergy to higher frequencies at the expense of the main 
super side-band mode; this process dominates for small 
to moderate initial steepness of the excited waves. 



•• Inverse cascade can be expected to be the prevailing 
process for high enough basic carrier wave steepness. 

•• Only a finite number of steps in both direct and 
inverse energy cascade leading to both breaking or stabi- 
lization regimes satisfy the BF-criterion. 

•• Fermi-Pasta-Ulam recurrent process may occur if 
at some cascade step modes are generated which are in 
exact resonance, i.e. AA; = in (30) (this is only neces- 
sary condition). This general prediction is in accordance 
with experimental results where it was shown that the 
observed recurrence is very close to the NLS or three- 
wave system of the form (18) and "neither a reduction in 
the number of waves per group nor downshifting of the 
spectral energy" is observed in the experiments of Tulin 
and Waseda ([5], p. 210) conducted in the laboratory 
tank of the sizes 50 x 4.2 x 2.1 m. 

The evolution of wave trains on an effectively much 
longer fetch than previous studies (330 x 5 x 5 m) was ex- 
perimentally investigated in [3, 4] with following results. 
A very long scaled wave modulation with several modula- 
tion loops demonstrates re-stabilization processes: finally 
the periodic modulations of wave train are observed at 
the latest stage with the most energetic lower side-band 
wave. This means the final termination of the wave cas- 
cade mechanism and stabilization of the system with es- 
sentially large steepness in accordance with our model. 

For small initial steepness of the wave packet £ ~ 0.1 
Benjamin-Feir instability leads to growing of the main 
pair side-band modes. Asymmetrical growth of side- 
bands with prevailing of the lower side-band mode and 
discretized energy cascade to the higher frequencies was 
also observed in the experiments [3, 4]. Near recurrence 
FPU phenomenon was clearly seen at the latest stages of 
the wave propagation: most part of wave energy revert 
back to the carrier frequency mode. For large enough 
initial steepness of the wave packet e ~ 0.15 to 0.25 wave 
steepness during propagation leads to wave breaking. Pe- 
riodic modulation and demodulation of wave trains are 
found at post-breaking stage, in which the energy of wave 
train transfers between the carrier wave and a pair of 
sidebands. Only partial FPU phenomenon was observed 
in this case - essential part of energy is lost due to break- 
ing. 

Necessary and sufficient conditions for manifestation 
of this phenomenon depend explicitly on the relation be- 
tween wavelength and initial phase of the carrier wave 
and the aspect ratio of the laboratory tank (see (17) 
and remark afterward; more details can be found in [6]). 
Deduction of their explicit form is outside the scope of 
this Letter where we concentrated on the study of non- 
periodic energy transfer due to the novel mechanism - 
dynamic cascade given by (24). 

• Our model of BF-instability can be further refined 
in many aspects. For instance, dissipation can easily be 
included at each step by changing of C„ to p„C„ in (31), 
with some constants < p„ < 1. Accordingly, the 
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form of the increments Idir.+m hm.-n will be modified 
and the number of cascade's steps will be reduced. De- 
tailed study of the effect of dissipation can be found in 
[11]. 

• This model is quite general and describes basic en- 
ergy cascade in an arbitrary 4-wave system with narrow 
initial excitation though the estimate (25) and definition 
of the small parameter e might change, depending on 
the specifics of the wave system. In particular, the 
model can be used directly as a basic description of 
wind generated instabilities of surface water waves from 
which extreme (or rogue) waves originate, e.g. see [12] 
for numerical simulations and [13, 14] for laboratory 
study. Theoretical study of the role of Benjamin-Feir 
instability in formation of extreme waves can be found 
e.g. in [15]. 
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